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A dilute two-dimensional (2D) Bose gas at zero temperature is studied by the method developed 
earlier by the authors. Low density expansions are derived for the chemical potential, ground state 
energy, kinetic and interaction energies. The expansion parameter is found to be a dimensionless 
in-medium scattering amplitude u obeying the equation 1/u + Inw — — ln(na^7r) — 27, where na'^ 
and 7 are the gas parameter and the Euler constant, respectively. It is shown that the ground state 
energy is mostly kinetic in the low density limit; this result does not depend on a specific form of 
the pairwise interaction potential, contrary to 3D case. A new form of 2D Gross-Pitaevskii equation 
is proposed within our scheme. 
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Theoretical investigation of 2D Bose gas is of inter- 
est not only in itself but also from the point of view of 
its experimental applications. Indeed, the experimental 
observation of the Bose-Einstein quasicondensate in the 
Bose gas of hydrogen atoms was claimed . On the other 
hand, the discovery of the Bose-Einstein condensation in 
magnetically trapped alkali-metal atoms |^ stimulated 
rapid progress in optical cooling and trapping of atoms. 
This progress gives us hope that an experimental obser- 
vation of 2D quasicondensation in the trapped atoms is 
a matter of near future (see some experiments along this 
line in Ref. |^). Theoretically, the Bose-Einstein con- 
densation is associated with the off-diagonal long-range 
order, i.e., the non-zero asymptotic at r = |ri — ^ 00 
(from the physical point of view at r ^ 1/ y/n) for the 
one-body density matrix 

(V'nri)^(r2)>-(^t(r,))(^(r2))^0. (1) 

Here tp''{r) and tp{r) are the Bose field operators, (• • •) 
stands for the statistical average, and (V'(r)) — 4>M is the 
order parameter. As was shown by Hohenberg H] from 
the Bogoliubov "l/g^" theorem in 2D case there is no 
off-diagonal long-range order at finite temperatures due 
to the temperature long-range fluctuations of the phase, 
and the limit (|^) is equal to zero. In spite of this fact, a 
phase transition is possible to a superfluid state at suffi- 
ciently low temperature Tc At this temperature, the 
asymptotic behaviour of the one-body density matrix at 
r ^ 00 is changed from an exponential decay (above Tc) 
to a power decay (below Tc) with respect to r. Thus one 
can speak about the phenomenon of the quasicondensa- 
tion in two dimensions. For zero temperature, the limit 
(^ differs from zero, and, hence, there exists the true 
Bose-Einstein condensate. Assuming that the condensate 
does exist at T = 0, in this paper we consider low-density 
expansions for 2D homogeneous Bose gas with respect to 
the gas parameter na^, where n = N/S is the density 
(a number of particles per unit area), and a stands for 
the two-dimensional scattering length (see Appendix ^). 



Note that the density expansion for the chemical poten- 
tial is intimately related to a form of the Gross-Pitaevskii 
equation, a powerful tool for investigating a dilute inho- 
mogeneous system of Bose particles . As the chemical 
potential is a continuous function of temperature, one 
can expect that our results, obtained for zero tempera- 
ture, are valid also for finite temperatures T <C Tc. Below 
we consider the case of zero temperature only. 

The leading term of the energy expansion in na^ for 
a gas of hard discs was first obtained by Schick ||], 
who made use of the Beliaev method 0, developed for 
3D Bose gas. Recently the Schick asymptotic formula 
has been proved rigorously To the best of our 

knowledge, there is only one paper, by Hines et al. [ pdl , 
where the next-to-leading terms were evaluated. How- 
ever, the authors of Ref. [pi]] employed the first-order 
Beliaev approximation for obtaining the next-to-leading 
terms, while in three dimensions the second-order Beliaev 
approximation is needed for the same purpose. Follow- 
ing Schick, they ignored the imaginary part of the Beliaev 
equation for the chemical potential (compare Eq. (4.3) in 
Ref. H with Eq. (1) in Ref. (n)). It is the unphysical 
imaginary correction to the chemical potential that de- 
termines the range of validity of the first-order approx- 
imation itself. One can easily demonstrate that in 2D 
case the correction is of the order of ri/ln^(na^); there- 
fore, that method can not yield correct terms of the 
expansion in this order and higher. At the same time, 
our method successfully reproduces in 3D case the 
famous next-to-leading term for the chemical potential, 
which ensures that our results are valid also in two di- 
mensions. Thus it becomes clear why only the first two 
terms involved in our expansion [see Eqs. (|l8|), ( p^ ) be- 
low] coincide with those of the corresponding expansion 
from Ref. ]|ll|: the correct coefficient even for the third 
term is beyond the approximation used in Ref. pj] ] . 

In this paper we adopt the method developed and 
described in detail in our previous publications [ p^ . 
Only basic notations and some important points are dis- 
cussed here. For a homogeneous system, the one-body 
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density matrix (■0'l'(ri)?/;(r2)) depends on r = ri — r2, 
and, hence, its eigenfunctions and eigenvalues are the 
plane waves exp(ip • v)/\/S and the occupation numbers 
Up = (a|jap), respectively. The Bose-Einstein conden- 
sate corresponds to the macroscopic occupation number 
Nq, and the order parameter is (^/'(r)) = {ao)/VS = 
^Jn^e^^ (uq stands for the density of the condensate). 
In turn, the eigenfunctions of the two-body density ma- 
trix (?/''l'(ri)-(/''l'(r2)'i/'(r2)-0(r'j)) can be naturally classified 
as follows. The maximum eigenvalue No{No — 1) ~ Nq 
corresponds to the state of two particles in the conden- 
sate; its eigenfunction (p{r) / S can be interpreted as a pair 
wave function in medium of the condensate-condensate 
type. The other macroscopic eigenvalues iN^Uq cor- 
respond to the two-body states with one particle in 
the condensate and another one beyond the conden- 
sate; its eigenfunctions (pq/2(r) exp[iq • (ri -I- V2)/2]/S 
are of the condensate-supracondensate type, where Sq 
is the total momentum of the pair of bosons. The 
residuary non-macroscopic eigenvalues are related to the 
supracondensate-supracondensate pairs and to bound 
ones provided the latter exist. The functions ip{r) and 
¥'q/2 (r) can be chosen as real quantities given by {p ^ 0) 

ip{r) = 1 -f ipp{v) = V2cos(p • r) + Vp(r) (2) 

with the boundary conditions ili{r), i/'p(i") — > at r — > 
00. The Fourier transforms of the scattering parts can 
be expressed in terms of the Bose operators: 
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(3) 



With the help of the in-medium scattering ampli- 
tudes U{k) = / d^r (y5(r)y(r) exp[— ik • r] and ?7p(k) = 
J d riy9p(r)F(r) exp[—zk-r], the chemical potential reads 



^l = noU(Q) + V2 



%c^q/2(q/2). 



(4) 



Here we introduce a pairwise interaction potential Vir). 
It should be emphasized that the formulas (||) and 
derived within the Bogoliubov principle of the corre- 
lation weakening are exact. For a Bose gas, a system 
with a small condensate depletion (n — no)/n ^ 1, the 
pair distribution function is expressed as 



V n / n J yl-Ky n ^' 



(5) 



where the contribution of the supracondensate-supracon- 
densate pair wave functions can be neglected. Another 
restriction for this representation is the assumption that 
there are no bound pair states In order to fulfil the 
latter condition, it is sufficient to require V{r) > 0, and, 
as usually, V{r) for r — > 00. In the framework of our 
scheme, the following equations are valid at sufficiently 
low densities fl^ : 



Tk + nU{k) 



f 2nTkU{k) 
U{k) 



- 1 



2nTkU{ky 



(6) 
(7) 



with Tfe = h'^k'^/{2m). Equations (|) and (|) look hke 
those of the modified Bogoliubov model where the "bare" 
pairwise potential V{k) is replaced by the effective one 
U (k) that is determined from the two-body Schrodinger 
equation. However, in our method there exists a key 
difference, which is of particular importance in two di- 
mensions: Eq. is a self- consistent equation for the 
in-medium scattering amplitude U{k). Indeed, using the 
definition of the latter, Eq. can be represented in the 
Lippmann-Schwinger form 



U{k) = V{k) - - 



1 f A\ l/(|k - q|)C/(g) 



(27r) 



T2 + 2nTqU{q) 



(8) 



Besides, one can make use of the limiting relation [see 
Eq. (|2|)] limp_,o ((5p (r) = ^/2ip{r), which allows one to 
simplify Eqs. (||) and (|) 



= n[/(0)[l 

5W = <^'W[n 



- {n - 
2(n 



m)/n- 

- no)/n- 



(9) 
(10) 



Thus, our scheme is reduced to the following. First, one 
should solve Eq. (||) and find U (0) [and (p{r)] as a func- 
tion of the density at n ^ 0. Second, the condensate 
depletion (n — no)/n should be determined from Eq. (|^). 
Third, employing these results, one should obtain the 
density expansion for the chemical potential and the 
short-range behaviour for the pair distribution function 
(p^. Note that in three dimensions this scheme is in 
excellent agreement with the data of Monte-Carlo calcu- 
lations for hard spheres (see the last paper in Ref. |]l2|). 

In order to solve Eq. (||) at n — > 0, we employ the 
procedure of linearization, which is similar to that in 3D 
case [E2I, and rewrite this equation in the form 



U{k) = V{k) - ip.P. 



d\ F(|k-q|)C/(g) / 



{2n) 



T —T 



2' (11) 



where P.P. denotes the Cauchy principal value, and / is 
given by 



/ = P.P. 



i2nr 



V{\k-ci\)U{q) n|k-q|)C/(g) 



r2 + 2nTqU{q) 



T —T 



Here we introduce the auxiliary quantity qq = 
c^j2mnU {0)/fi, where c stands for an arbitrary dimen- 
sionless constant. Performing the "scaling" substitution 



ci'y/2mnU{0)/h 



(12) 



in the integral [U{0) is assumed to depend on n in such a 
manner that nU (0) — > when n 0] and, then, taking 
the zero-density limit in the integrand, for n — *■ we find 
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Then, with the help of the Fourier transformation, 
Eq. (|ll|) reads 



Thus, with the help of Eqs. (|T|), (|T|) and (|ol), one can 
rewrite Eqs. and ( pO| ) as 



1- A 



4^2 



d\'V{r')^{r')Y^{q^\v-v'\). 



fi = {ATrh^n/rnjull + u + 



(21) 
(22) 



Here the Fourier representation is used Yo{pr) = 
4/(27r)2p.P./ d^q exp[iq • r]/(p2 _ g2) ^j^g cylindri- 
cal Bessel function of the second kind. Since only the 
asymptotic of ip{r) at n ^ is of interest, the linear 
integral equation for ip{r) can be written as 



(fir) 



1- A 



Note that Eqs. (|l^) and ( ^2|) are the short-range approx- 
imation valid at r < 1/ ^/n. For this reason, the bound- 
ary condition (||) is not fulfilled for ip{r) in Eq. ([T^). 
In order to obtain the energy per particle e, we rep- 
resent it in the form e — {2TTlT?n/'m)f{u) with an un- 
known function f{u). From Eq. ([l^) it follows that 
ndu/dn — m^/(1 — u), which, together with the thermo- 
dynamic relation fi — d{en)/dn, yields the differential 
dVy(r')^(r') ln(qo|r - r'\e^/2), equation u2d//du + 2(1 - u)f = 2u{l - u^). From this 

equation and the initial condition f{u = 0) we derive 



where the asymptotic Yo{z) = 2 ln(2:e'^/2)/7r -I- 0(2^ In z) 
for z ^ is used. Here 7 ~ 0.5772 stands for the Euler 
constant, and 0{x) denotes terms of the order of x or 
even higher. It is seen from the resulted equation that, 
first, ip{r) obeys the Schrodinger equation (Al), and, sec- 
ond, its asymptotic for r — > cx) is 



ip{r) 1 - A + ln(rgoe'^/2)m(7(0)/(27r;i2). 



(14) 



which differs from that of Eq. ( |A2| ) only by the multipli- 
cation factor m[/(0)/(27r/i2). Compari ng E q. (|^ with 
Eq. ([A^) yields due to linearity of Eq. (Al) 



ip{r) =2uip^°^{r), (15) 
-Ina = 2TTh'^{l - A)/[mC/(0)] -I- ln(goe''/2), (16) 

where we introduce the parameter u by the definition 



U{0) = / d'^rip{r)V{r) = {4Trh^/m)u. 



(17) 



With the help of Eq. (|l^ and the definition of qq (see 
above), Eq. ( |l6|) can be rewritten as 



(5(1 + ulnu) 



-l/[ln(na27r) + 27]. (18) 



Note that (5 ^ at n — > 0. As expected, the arbitrary 
constant c is cancelled and not involved in final Eq. ( 1^ ) 
and, hence, in the formula ( p^ for if{r). Equation (1^) 
has no solution for u when 5 > 1 and has two positive 
ones when S < 1 (i.e., when na^ < 0.0369 . . .). The solu- 
tion with a greater value of u should be ignored because 
of its unphysical behaviour [u ~ l/(na^ at n — > 0]. An 
expansion for u is obtained from Eq. ( |l8| ) by iterations 



u = 5 + S-'lnS + S-'lnU + S^\n5 + OiS^ln^ S). (19) 

Using Eq. (|^), one calculates the condensate depletion 
by means of the substitution ( p^ upon integrating 



no 



(20) 



£ ^ {2TTn^n/m)[u + u^/2 + 0{u^)]. 



(23) 



With Eq. (^9|), e can be expanded in the parameter S 

e— {2nh^n/m.) 
x[S + 6^ \n5 + S^/2 + 5^ In^ S + 2S^ \nS + 0{S^)]. (24) 

One can see that all the low density expansions are series 
in the dimensionless in-medium scattering amplitude u, 
which depends ultimately on the density via Eq. ( p^ ) 
[and, hence, ([l9|)]; therefore, u can be considered as a 
parameter of low density expansions in two dimensions. 
The interaction energy per particle is exactly related to 
the pair distribution function ( ^2|) 



£int 



d\V{r)g{r) = ^^^^[u' + 2u' + ■■■], (25) 



where we put by definition 



(26) 



In the latter equation we employ the theorem (A4) with 
the coupling constant A [V{r) AF(r), and A = 1 
in final formulas!. Since ei^it can be directly evaluated 
via the relation (Eq), our approach takes accurately into 
account the short-range particle correlations jl^]. Note 
that Eq. (^5|) can also be obtained from the Hellmann- 
Feynman theorem Si^t = Xde/dX with Eqs. (|l^), (|2^), 
and (p6|). Moreover, Eq. ( |2^ ) can be derived in the same 
manner varying the energy: g{r) = [2/n)5e/5V{r). For 
the kinetic energy per particle Ckin — {J2 Pf) / i'^^-^) '^^ 
have 



2Trh'^n r 



^int 



(27) 



It is seen from Eq. (^3|), ( p5| ) and (|27| ) that in the lead- 
ing order, proportional to nu, the total energy is purely 
kinetic. Thus, whatever a particular shape of the po- 
tential V{r), at sufficiently small densities the energy 
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becomes mostly kinetic. By contrast, in three dimen- 
sions £kin ~ 2iTh^bn I m and ej^t ~ 21111? {a — h)n/ra are of 
the same order, where a is the 3D scattering length and 
h = a — Xda/dX (except for hard-spheres when a = b, 
see Ref. |3|). 

Let us discuss the nature of the Schick approximation 



-2TTf?n/{m\n nc? ) , 



(28) 



which is Eq. (EJ) in the lowest order in njhinc?. As 
the energy (28) in this order is purely kinetic, it can- 
not in principle be represented as a sum of the interac- 
tion energies of two particles over all pairs of bosons by 
analogy with the weak-coupling 3D Bose gas. However, 
as in three dimensions, we can start from Eq. (^ and 
put /i ~ n J d^r V^(r)(p(r) with the in-medium pair wave 
function <y3(r). It is clear that 95 (r) for r < (here 
ro ^ is of the order of the correlation length) 

should be proportional to the wave function (^(°)(r) of 
the two-body problem (^): Lp{r) ~ C<f^°\r). The 
boundary condition (|^) can be fulfilled only due to in- 
medium effects for r > vq; therefore, we can approxi- 
mately put f{ro) ~ 1 and use for 3> a the asymptotic 
(^). This leads to C ~ -2/lnna2 and, by E^^ (^), 
yields /i ~ —Airh^n/ (mhina?) and, hence, Eq. (psj). The 
crucial difference in 3D case is the boundary condition 
for the two-body problem ip'^'^\r) ^ 1 — a/r instead of 
Eq. (|A2|). In this case the condition iy9(ro) ~ 1 leads 



to ip{r 




(r) in the leading order, which results in 



/i ~ rt / d^r V{r)(p'''^'> (r) = A-Kh^na/m. 

Now one easily writes 2D Gross-Pitaevskii functional 
for the energy using Eq. (^) in the leading order 



dVf^l^^ 



2m 



+ Vcxt(r) 



27rr 



(29) 



and 2D Gross-Pitaevskii equation 

ihd(l)/dt = SE/S(t>* 
= [-(ftV2m)V2 + Vc^t(.r)](j) + {'inf?u/m)\cb\^<p. (30) 

Here (j> = 0(r,t) = {ijj{r,t)) is the order parameter with 
the normalization N — j"d^r|0p, and u is given by 
Eq. ( p^ with n ~ \4>^ ■ Upon varying in Eq. (pO[), we ne- 
glect the variation of m, for \^\^8ul ^ 'u?\^<\). Note 
that the sum of the first and third terms in Eq. ( |29| ) 
corresponds to the kinetic ene rgy o f bosons according to 
Eq. (|2^). Equations (|2S| ) and ( pO|) are more exact than 
those of Ref. jl^ , based on the Schick formula (^8|) . 

In conclusion, the expansions have been derived for the 
condensate depletion (pO|), the chemical potential ([2l]), 
the pair distribution function ( p2|) for r < 1/ ^/n, the to- 
tal (p3|), interaction ( p5| ) and kinetic ( p7| ) energies. The 
energy expansion (^3|) leads to the Gross-Pitaevskii func- 
tional (E3) and equation (30). This work was supported 
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APPENDIX A: 

In this appendix, a useful variational theorem is proved 
for 2D scattering length. In this paper we deal only with 
the short-range potentials which go to zero for r — > 00 
as V{r) 1/r™ (m > 2), or even faster. The s-wave 
function corresponding to relative motion of two particles 
with p = obeys the two-body Schrodinger equation in 
the centre-of-mass system 

(°)(r)-|-y(r)(^(")(r) =0 (Al) 

(0)(r)| < 



with the following boundary conditions: first, \ip 
00 at r = 0, and, second, for r — > 00 

(/j(°)(r) ^ ln(r/a). 



(A2) 



Since the asymptotic is chosen to be real, the solution of 
Eq. (Al) is also real. The introduced positive quantity a 
is called 2D scattering length. Integrating Eq. ( [Al[ ) and 
keeping in mind Eq. (|A2|) yield 



(A3) 



Let us suppose th at V(r) is infinitesimally changed. 
Then, varying Eq. (Al), multiplying the obtained equa- 
tion by ip^^\r), and carrying out the integration, one 
arrives at the theorem 



2Trt? Sa 



= / dV[^(")(r)]25V^(r). 



(A4) 
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